function f = AR1_Adj(rho , T , order)

if order == 1
    f = rho + (1 + 3 * rho)/T;
elseif order == 2
    f = rho + (1 + 3 * rho)/T + 3 * (1 + 3 * rho)/T^2;
end